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Abstract 

Power-suppressed corrections coming from the end-point integration regions in 
the amplitude of the process 7*7 — > nn at large Q 2 and small squared center- 
of-mass energy W 2 are calculated in the QCD hard-scattering approach where 
the amplitudes factorize in a hard perturbatively calculable part and a general- 
ized distribution amplitude. The running coupling method and the technique of 
infrared renormalon calculus are applied to obtain Borel resummed expressions 
for the two main components of the process amplitude. Numerical estimates for 
these power corrections are presented. They are sizeable when Q 2 < 10 GeV 2 . 



1 Introduction 

Single meson and meson pair productions with a small invariant mass W in virtual photon- 
photon collisions are exclusive processes for which perturbative QCD (PQCD) [1] analysis 
was successfully applied when the virtuality Q 2 of one photon is high. These investigations 
are based on the perturbative QCD factorization theorems which allow one to compute the 
amplitude of the exclusive process as the convolution integral of the meson distribution 
amplitude (DA) or the two-meson generalized distribution amplitude (GDA) 2 4 and the 
hard-scattering amplitude of the underlying partonic subprocess. The meson DA's and 
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GDA's <p(x,fi 2 F ) and <f>(z, (, W 2 , fip) are non-perturbative objects and contain the long- 
distance mesonic binding and hadronization effects. The GDA's are related by crossing [3 
to generalized parton distributions (31 |3] . The theoretical and experimental investigations 
of the processes 7*7 — > MM open opportunities to obtain new, valuable information on 
the fragmentation of quarks or gluons into mesons |6 r These processes can be studied in 
ey and e + e~ collisions and first results have been published [7]. 

The perturbative QCD approach and the factorization theorems describe exclusive 
processes at asymptotically large values of the squared momentum transfer Q 2 . But 
in the present experimentally accessible energy regimes, power-suppressed corrections 
may play an important role. There are numerous sources of power corrections to the 
process 7*7 — > MM. For example, power corrections arise due to the intrinsic transverse 
momentum of partons retained in the corresponding subprocess hard-scattering amplitude 
and GDA's. They may be estimated along the lines presented in Ref. jEJ, where such 
corrections were calculated for the deeply virtual electroproduction of photon and mesons 
on the nucleon. A power correction to the process 7*7 — ► ttti has been estimated for 
the amplitude corresponding to scattering of two photons with equal helicities, with the 
help of the light-cone sum rules method [§] . The power-suppressed (twist-3) contribution 
due to the interaction of a longitudinally polarized virtual photon with the real one was 
analyzed within the Wandzura-Wilczek approximation in Ref. [TO] . 

In the present paper we compute a class of power corrections which originate from the 
end-point regions z — > 0, 1 in the integration of the PQCD factorization expression over 
the parton longitudinal momentum fraction z. We restrict ourselves to the two-pion final 
state and to the leading twist- 2 amplitudes. Generalization of our approach to encompass 
other two-meson final states is straightforward. 

It has been advocated JT] that, in order to reduce the higher-order corrections to a 
physical quantity and improve the convergence of the corresponding perturbation series, 
the renormalization scale, i.e. the argument of the QCD coupling in a Feynman diagram 
should be set equal to the virtual parton's squared four-momentum. In exclusive processes, 
the scale chosen this way inevitably depends on the longitudinal momentum fractions 
carried by the hadron constituents. In our case, the relevant scale is given by = Q 2 z 
or Q 2 z 0. But then the PQCD factorization formula diverges, since a s {Q 2 z) [a s (Q 2 z)} 
suffers from an end-point z — > [z — > 1] singularity. This problem may be solved by 
freezing the argument of the QCD coupling and performing all calculations with a s (Q 2 ) 
[or a s (Q 2 /2)]. In the running coupling (RC) method, one allows the argument of a s to run 
but removes divergences appearing in the perturbative expression with the help of a Borel 
transformation and a principal value prescription. It turns out that this procedure, used 
in conjunction with the infrared (IR) renormalon technique ^21 EH] allows one to obtain 
the Borel resummed expression for the process amplitude and estimate power corrections 
arising from the end-point integration regions. This method was used for other processes 
[HJ [13 EH H3 CHj and succesfully confronted to experimental data. 

This paper is organized as follows: in Sect. 2 we present kinematics, general expres- 
sions for the amplitude of the process and the two-pion GDA's. In Sect. 3 we outline 
the main points of the RC method and obtain the Borel resummed components of the 
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Figure 1: Schematical representation of the factorization theorem for the process 7*7 — > 
7T7T. The solid and dashed blobs denote the hard-scattering subprocess 7*7 — > gg (7*7 — ► 
and hadronization qq —>■ (gg — > 7T7t), respectively. The quarks (gluons) are depicted 
as dashed lines. The momenta of the initial and final particles are shown in the figure. 
The total momentum of the final state is P = p\ +p 2 and z is the longitudinal momentum 
fraction of the 2it system carried by the quark (gluon). 

amplitude. Section 4 contains results of our numerical calculations. Finally, we give our 
conclusions in Sect. 5. 

2 Amplitude of the process 7*7 — > tttt 

The process 



is schematically depicted in Fig.^ In addition to the four-momenta of the initial and final 
particles one introduces the total and relative momenta of the pion pair P — pi +P2, A = 
P2~P\- The momenta of the involved particles can be described in terms of two light-like 
vectors p, n which obey p ■ n = 1. The decomposition of the four- momenta of the initial 
and final states in terms of the vectors p and n reads 



7*(?)+7(?') 



7r(pi) + n(p 2 ), 



(1) 



Q 2 , 

—n, q = 



Q 2 + W 2 
2 



n, q 2 = -Q 2 , q' 2 = 0, 



W 2 - - w 2 

P=P + —n, P 2 = W 2 , A = ((-Op+((-0—n + k ± , (2) 



where the quantities 



C = ^— , C = i-C 



P ■ n ' P • n 

describe the distribution of the longitudinal momentum between two pions. The vectors 
p and n can also be employed to define the metric tensor in the transverse space: 

{-g^) T = -g^ + p*W + p v n^. (3) 

In the hard photoproduction regime Q 2 3> W 2 , A 2 the amplitude of the process (JTJ) has 
the form 



r*"(C, w 2 ) = i (-n T r (c, ^ 2 ) + i ^+^ r^c w^ 2 ) + ^^r 2 (c, w^ 2 ), (4) 

where k±k± is the traceless, symmetric tensor product of the relative transverse momen- 
tum of the pion pair 

k±k ± = kj_k± — - (—g^ u ) T k\. 

In (jlj) To((,W 2 ) is the amplitude corresponding to the scattering of two photons with 
equal helicities, Ti(C, W 2 ) denotes the amplitude with L z = ±1, whereas T 2 (£, W 2 ) arises 
from the subprocess with opposite helicity photons. In fact, in the collinear approxima- 
tion, conservation of the angular momentum along the collision axis leads to the helicity 
conservation h* — h = h\ + h2, where h*, h are the helicities of the virtual and real photons 
and hi, h 2 denote the helicities of the produced quarks or gluons. When photons pro- 
duce a quark-antiquark pair, at the leading twist-2 level only the subprocess with L z = 
contributes to the amplitude. The subprocess with L z = ±1 is twist-3 and determines 
the contribution T%((, W 2 ) appearing due to the interaction of a longitudinally polarized 
virtual photon with the real one ^U] • I n the case where photons create a gluon pair, both 
the subprocesses with L z = and L z = ±2 contribute at the twist-2 level to (J3J). 

We now focus on the amplitudes W 2 ), i = 0, 2 which do not vanish at leading 
twist. They can be written as convolution integrals of hard-scattering coefficient functions 
C(z,fi 2 F ) and two-pion GDA's <&(z, (, W 2 , /j 2 f ), i.e 

To{(, W 2 ) =E e ? f dzC q {z, » F )%(z, C, W 2 , /4)- 

j 

£ e\ f dzC g (z, » 2 F )<S> g (z, C, W 2 , /j, f ) (5) 

and 
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Figure 2: Sample Feynman diagrams of the hard-scattering subprocesses 7*7 — > gg (at 
leading order- a, at next-to- leading order -b) and 7*7 — ► gg (c). 



T 2 (C,^ 2 ) = J^dzCjiz^D^^CW 2 ,^, (6) 

where /i^ is the factorization scale at which the hard and soft parts of the reaction are 
defined. Here the coefficient function C q (z,fj, F ) is calculated with NLO accuracy using 
the subprocess 7*7 — > qq (Figs|2K, b) with L z = 0, whereas the functions C g (z,fi 2 F ) and 
Cj(z,fj, 2 F ) correspond to the subprocess 7*7 — > gg (Fig. with L z = and L z = ±2, 
respectively, and contribute only at the next-to-leading order of PQCD due to quark-box 
diagrams. 

The amplitudes Ti((,W 2 ) do not depend on the renormalization and factorization 
schemes and scales employed for their calculation. But at any finite order of QCD pertur- 
bation theory, due to truncation of the corresponding perturbation series, the coefficient 
functions depend on both the factorization \j? F and renormalization scales. An optimal 
choice for these scales is always required to minimize higher order corrections. The fac- 
torization scale n 2 F in exclusive processes is traditionally set equal to the hard momentum 
transfer Q 2 , and we shall follow this prescription. 

The functions C(z) with NLO accuracy are given by the following expressions JH] 



C q {z)=C° q (z) + ^-C 1 q (z), (7) 



C g (z) = °^kc]{z), Cj(z) = ^C T (z), 



and 



r ln 2 z In 2 z In 2 z In 2 ~z In z lnz 9 9 
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Co ( z ) = o_ 9 z 2 In 2 z + z 2 \n 2 z + 2zz In zz — Az In z — Az In z 

y z i z i l 

C\{z) = 4, (9) 

zz 

with CV = 4/3 being the color factor. 

Choosing the renormalization scale yr R is more subtle and we follow the prescription 
[TT] that it is equal to the square of the momentum transfer carrying by a virtual parton in 
each leading order Feynman diagram of the underlying hard-scattering subprocess. Here, 
these scales are determined by the leading order diagrams of the subprocess 7*7 — > qq 
and are given by the virtualities of the off-shell fermion, which are equal to Q 2 z or to Q 2 z 
depending on the diagram. In the present paper we adopt the symmetrized RC method, 
where a s (Q 2 z) and a s (Q 2 z) are replaced by 

a s {Q z), a s {Q z) => . (10) 

The reasons which led to introduction of this version of the RC method and further details 
have been presented in [T7| ITHj. 

The next component in the factorization formulas (JSJ), © is the generalized distribu- 
tion amplitudes of the 2tt system. At present, little is known about these GDA's, but 
using constraints originating from crossing symmetry and soft pion theorems, as well as 
the evolution equation for GDA's, we can model them. Further information on their 
analytic form should be extracted from analysis of experimental data and corresponding 
theoretical predictions and/or obtained, as in the case of the usual DA's of mesons, via 
QCD non-perturbative methods. 

The simplest 2n GDA's obtained using the requirements described above are 

$ q {z, C, W 2 , /4) = 20zz{z - z)—M q (^ F )A(C, W 2 ), 

%(z,(,W 2 ,li 2 f ) = 10z 2 z 2 M a (n 2 F )A((,W 2 ), (11) 

where the model-dependent (201 function A((, W 2 ) is z and Q 2 independent. Our results 
will not depend on the choice of this function. 

The GDA's represented by the formula (jllj) depend on the momentum fractions carried 
by quarks M q (n 2 p) and gluons M g (fip) in the pion 

MM) = MTU + R{tiW F )}, R(fi) = ^ q(l4) M:r 



asy 



with 



L(i4 



« s (/4) 



Mi 



7+ = l(n f + 4CV), M g ^ 2 F ) = 1 - MM), (12) 



a s (Mo) 

the asymptotic values of which are determined by the expressions 
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lyjasy _ n f M asy = ^ F 

q " n f + 4C F ' 9 ~n f + AC F 



For the helicity-two GDA <t^(z, C, W 2 } up) we take 



*T(z, C, w\ nl) = d^ 2 f )z 2 z 2 a t u, w 2 ) 



(13) 



(14) 



with 



£>G4) = 30D T (/ig) 



a s (// 



«s(>o) 



00 



7 



TG 



7+- n/ . 



(15) 



In (JEJ), (fTH|l /1q and DT(/Jq) are the normalization point and constant, respectively. 



3 Borel resummed amplitudes 

Computation of the amplitudes Ti(Q 2 , (, W 2 ) implies integrations over z. Having in- 
serted the explicit expressions of the hard-scattering coefficient functions and the two-pion 
GDA's into (jSJ and (jOJ we encounter divergences, arising from the singularities of the cou- 
pling constant a s (Q 2 z) and a s (Q 2 z) in the limits z — ► 0, 1. The RC method proposes a 
way to cure these divergences. 

To this end we express the running coupling a s (Q 2 z) in terms of a s (Q 2 ) [a similar 
argument holds also for a s (Q 2 z)}. This is achieved by applying the renormalization-group 
equation to a s (Q 2 z) [21]. We get 



« s (Q 2 )/3iln[l + ^]l 



2tt/? 1 + ^ 



Here a s (Q 2 ) is the one-loop QCD coupling constant, t = 4n / j3 a s (Q 2 ) = ln(Q 2 /A 2 ) and 
/3q, (3\ are the QCD beta function one- and two-loop coefficients, respectively, 

2 19 

A) = 11 - gTl/, A = 51 - yTl/, 

and ?7./ is the number of active quark flavors. Equation (|TB|) expresses a s (Q 2 z) in terms 
of a s (Q 2 ) with an a 2 (Q 2 ) order accuracy. 

Inserting (jlfij) into the amplitudes we obtain integrals, which can be regularized and 
calculated using the method described in [T3j. The amplitudes Ti(Q 2 , £, W 2 ) are then writ- 
ten as perturbative series in a s (Q 2 ) with factorially growing coefficients C n ~ (n — 1)!. 
Their resummation is performed by employing the Borel integral technique [22] . Namely, 
one has to determine the Borel transforms B[Ti(Q 2 , (, W 2 )](u) of the corresponding se- 
ries and in order to find the resummed expression for the amplitudes, has to invert 
B[Ti(Q 2 X,W 2 )](u) to get 
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res r°° 




~ P.V. / duexp 
J Jo 





B 



T t (Q 2 ,C,W 



+ 



amb 



(17) 



The Borel transforms B [T^Q 2 , (, W 2 )} (u) contain poles {u } located at the positive 
u axis of the Borel plane, which are exactly the IR renormalon poles. Therefore, the 
inverse Borel transformation can be computed after regularization of these pole singular- 
ities, which is achieved through a principal value prescription. But the principal value 
prescription itself generates higher twist (HT) ambiguities (uncertainties), which in the 
right-hand side of (|T7|) are denoted by [T^Q 2 , (, W 2 )] amb . They are determined by the 
residues of the Borel transforms at the IR renormalon poles q G {uq} and depend also on 
unknown numerical coefficients {iV^} 



T t (Q 2 ,(,W< 



amb 



?6{uo} 



N, 



Q 2q 



The ambiguity (fTKj) can be used to estimate higher twist corrections to the amplitudes 
stemming from other sources [for example, from the 27r higher twist GDA's]. 

A useful way to avoid the intermediate operations and obtain directly the Borel re- 
summed expressions is to introduce the inverse Laplace transformations of the functions 
in (021), i.e. 



1 



1 



(t + x) 1 



r> Jo 



POD 

/ duexpl— u(t + xMu"' 1 , Rev > 0, 
Jo 



(19) 



and 



ln[t + x] 



/ duexpl— u(t + x)](l — 7^ — \nu)u, (20) 
(t + x) 2 Jo 

where T(i>) is the Gamma function, 7^ ~ 0.577216 is the Euler constant and x = \nz 
[x = \nz in the case a s (Q 2 ~z)\. Then, the QCD coupling a s (Q 2 z) may be written as 16 j 



477- poo 

a s (Q 2 z) = — / due- ut R(u,t)z- u , 
Po Jo 



with 



R(u,t) 



Po 



In t — In u) 



(21) 



(22) 



The expression for the QCD running coupling (|2T| is obtained from (j!6)l and is suited 
to account for the end-point effects. It differs from that introduced to perform the re- 
summation of diagrams with quark vacuum insertions ( "bubble chains" ) into a gluon line 
[121 US]- In exclusive processes both these sources lead to power corrections. As noted 
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above, in the present work we consider contributions to the process amplitudes arising 
only from the end-point regions. 

Using (|14|) . (J21j) and performing the integration over z we get 2 



\UQ\ C, W*)r = 20A((, W 2 ) £ e 2 q { f 1 + ^ r due^R(u, t)Q(u 

fdue- ut R(u,t)G(u)\, 



and 



with 



2/3 Jo 
D{Q 



(23) 



[T 2 (Q 2 ,C,W 2 )] res = 2A T g ^W 2 )^f>- ndue- ut R(u,t)B(2-u,2), (24) 

Po Jo 



Q{u) 



! + £,B(2,f3)\^ u - ^B(l - u,/3)\ 2 - J^B(1,/?)| 2 _ U + 



03 



dp 2 



W^B(2-u, 

^B(l-u,f3)\ 3 + ^B(l,(3)\ 3 . u -^B(2-u,(3)\ 2 -^B(2^)\ 2 _ u 
3 d_ B (l - u, p)\ 3 + 3±B(1, /3)| 3 _ - 3^5(2 - u 



3^5(2, (3)\ 2 _ u - 



95(3 - u, 1) - 95(1 - u, 3) + 185(2 - u, 2), 



(25) 



and 



G(«) = ^B(l - u, [3)\ 3 + £B(1, /3)| 3 _ w + 2^5(2 - u 



d(3 2 



'df3 J 



+ 



2^5(2, /3)| 2 _„ - 4^5(1 - it, /3) | a - 4^5(1, (3) | 2 _„, 



(26) 



where B(x,y) is the Beta function B(x,y) = T(x)T(y)/T(x + y). 

In order to proceed one has to reveal the IR renormalon poles in the resummed 
expressions. The analysis of the pole structure of Q{u) and G{u) is straightforward. 
The result is that the function Q(u) contains a finite number of triple poles located at 
u = 1, 2, 3, an infinite number of double poles at the points u = 2, 3, 4 . . . oo and 
single ones at the points u = 1, 2, 3, 4 . . . oo. For the function G{u) we get: triple pole 
with location at uo = 3, infinite number of double (wo = 2, 3, 4 . . . oo) and single poles 
(wo = 1, 2, 3, 4 ...oo). The amplitude [T 2 (Q 2 , (, W 2 )] res possesses only single poles at 
uq = 2, 3. In other words, by employing (|T6*j) we have transformed the end-point diver- 
gences in (JSJ) and © into the IR renormalon pole divergences in (|23|) and (J24)) . One 
can see that the integrals in these expressions are the inverse Borel transformations (|17|) . 
where the Borel transforms B q ^[T (Q 2 , £, W /2 )]('u) of the quark and gluon components 
of the amplitude [T (Q 2 , (, W 2 )] res (in the quark case the NLO part) and that of the 
amplitude [T 2 (Q 2 , (, W 2 )] res up to constant factors are defined as 

2 Below, for brevity, we do not write down explicitly the higher twist ambiguities in the resummed 
expressions. 
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B q[g] [T Q (Q 2 , C, W 2 )}{u) ~ t)Q{u) [-G(u)\ , B[T 2 {Q\ (, W 2 )]{u) ~ *)fl(2-u, 2). 



After removing IR renormalon divergences from (f2*3*J) and by means of the princi- 
pal value prescription, they determine the resummed amplitudes [Ti(Q 2 , (, W 2 )] res . The 
final expressions [Ti(Q 2 , (, W 2 )] res contain power-suppressed corrections ~ 1/Q 2n , n = 
1, 2, 3 ... to the amplitudes fSJ El EH] and are the main results of the present work. 

Let us now check the asymptotic limit of the resummed amplitudes. In the asymptotic 
limit Q 2 — ► oo , GDA's ® q (z, (, W 2 , Q 2 ) and $ g (z, (, W 2 , Q 2 ) evolve to their asymptotic 
forms obtainable from (jllj) by means of the replacements M q (Q 2 ) -> M™ y and M g (Q 2 ) -> 
Mg Sy . We need also to take into account that in this limit the sub leading term in the 
expansion of a s {Q 2 z) through a s (Q 2 ) should be neglected, i.e. 

/•oo roo 

/ due- ut R(u,t) / dwe""*. (27) 
Then the amplitude [T (Q 2 , (, W 2 )] res takes the following form 



Er oW ., C, «^)r - E«J {l + f f *^ - **0l} . (28) 

The asymptotic limit of the integrals can be computed using techniques, described in a 
detailed form in ^7J EI] ■ After some manipulations, one gets for the asymptotic limit of 
the amplitude [T (Q 2 , C, W 2 )] res 

T (O 2 C W 2 )] ^ - 2 ° A ^ W2) V e 2 1 1 87 r as{Q2) \ (20) 

This expression coincides with the corresponding result from ^J] 3 and can be read- 
ily obtained within the standard approach employing the 2ir asymptotic GDA's. The 
asymptotic limit of the amplitude [T 2 (Q 2 , C, V4 /2 )] res , due to the factor D(Q 2 ), is equal to 
zero. 

This analysis of the asymptotic limit of the Borel resummed amplitudes shows the 
internal consistency of the RC method itself. 



4 Numerical results 

Let us now present numerical estimates of the power corrections to the amplitudes. The 
resummed amplitude [T (Q 2 , (, W 2 )] res , which contains both the hard perturbative com- 
ponent and power corrections, can be rewritten in the form 

iUQ 2 ,(,W 2 )Y es = T LO (Q 2 ,( , W 2 ) + T NLO (Q 2 , C, W 2 ) + T PC (Q 2 , (, W 2 ), (30) 

3 Our definition of the function A((, W 2 ) differs by a factor —1/6 from that of |19| . 
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Figure 3: The ratios R 2 (Q 2 ), R(Q 2 ) and l+Ri(Q 2 ) as functions of Q 2 . The solid (dashed) 
curve corresponds to the parameter M q (l GeV 2 ) = 0.6 [M g (l GeV 2 ) = 0.5]. 



where the first two terms in the RHS of (jBllj) are the LO and NLO parts of the amplitude, 
whereas Tq c (Q 2 , (, W 2 ) denotes the power corrections to it. The latter is given by the 
expression 



res 



Trr c (Q 2 X,W 2 ) = [T (Q 2 ,(,W 2 )\ nlo -T n ^(Q 2 X,W 2 ). (31) 

For our purposes it is convenient to normalize the expressions (I30J) and ()31|) in terms of 
Tq°(Q 2 , (, W 2 ), which results in ratios independent on the function A((, W 2 ): 

R{Q 2 ) = 1 + R 1 (Q 2 )+R 2 (Q 2 ). (32) 

Here 

R(Q 2 ) T « LO W^ W2 ) T^(Q 2 ,CW 2 ) 

Ri{Q )- t lo (q2; c> W2) > «»W ) ~ T L 0( g 2j C W 2 } ■ W 

In our calculations we use the following values of the parameters A and fir, 

A 4 = 0.2 GeV, fi 2 = 1 GeV 2 . (34) 

To clarify the sensitivity of the predictions to the parameter M g (/ig) we shall take the 
two plausible values M q (l GeV 2 ) = 0.5 and M q (l GeV 2 ) = 0.6. We use the two-loop 
approximation for the QCD coupling a s (Q 2 ). 

The amplitude [Tq(Q 2 , (, W 2 )] res contains an infinite number of IR renormalon poles. 
We truncate the corresponding series at some n max = 50 which is amply sufficient [TTjlTH]. 
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Figure 4: The ratio R%(Q 2 ) vs Q 2 . The solid line corresponds to R^{Q 2 ) without the HT 
ambiguities. The broken lines are obtained by taking into account the HT ambiguities 
([IBD . For the dashed line: N 2 = N 3 = 1; for the dot-dashed line: N 2 = N 3 = -1. 

In Fig. |3]we show R 2 (Q 2 ) as a function of Q 2 - The power corrections amount to some 
50-60 per cent of the corresponding leading order contribution at Q 2 = 1 GeV 2 . They 
are not completely negligible also at Q 2 = 10 GeV 2 reaching around 15% of the LO term. 
One observes that in the region 1 GeV 2 < Q 2 < 4 GeV 2 the function R 2 {Q 2 ) is more 
sensitive to the chosen value of the parameter M q than in the domain Q 2 ~ 10 GeV 2 . 

In Fig.EJthe ratios R(Q 2 ), l + Ri(Q 2 ) are also shown. As is seen the power corrections 
significantly reduce the amplitude T (Q 2 X,W 2 ) and this effect depends on the 2n GDA 
used in calculations. Thus, at Q 2 = 1 GeV 2 the resummed amplitude computed using 
the 2-7T GDA with the input parameter M q = 0.6 is approximately twice as large as the 
same amplitude found employing the GDA with M q = 0.5. At the higher values of Q 2 
this difference becomes more moderate ~ 1.1 at Q 2 = 10 GeV 2 . 

Another conclusion, which can be made after analysis of Fig. El is that the difference 
between the resummed [the ratio R(Q 2 )] and the standard predictions for the amplitude 
[the ratio 1 + Ri(Q 2 )} becomes smaller at higher values of the momentum transfer Q 2 . 
In fact, at Q 2 = 1 GeV 2 the resummed amplitude is equal to 0.41 of the standard result, 
whereas at 10 GeV 2 one gets 0.85 [for M q (l GeV 2 ) = 0.6]. 

For phenomenological applications it is useful to parametrize the ratio R 2 {Q 2 ) using 
the power-suppressed terms ~ 1/Q 2n , n = 1,2,3. Our fitting procedure leads to the 
following expressions 



R 2 (Q 2 )~ -1.709 + 



1.881 



0.7524 



M q (l GeV 2 ) = 0.5 



Q 2 



Q 4 
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-1.462 + 



1.515 0.533 



Q' 2 



Q 1 



M q {l GeV 2 ) 



0.6. 



The power corrections to the amplitude T 2 (Q 2 , (, W 2 ) are given by the formula 



T 2 PC (Q*,CW< 



The ratio 



T 2 (Q 2 ,C,W 2 . 



T 2 (Q 2 ,(,W 2 



(35) 



(36) 



T 2 (Q 2 ,C,^ 2 ) 

is shown in Fig. 0] It turns out that in this estimate of the power corrections to the 
amplitude T 2 (Q 2 , (, W 2 ) are large and may still amount to a 60 per cent increase of the 
amplitude at Q 2 ~ 10 GeV 2 . Such a large magnitude of the end-point effects can be traced 
back to the fact that T 2 (Q 2 , (, W 2 ) begins at 0(as(Q 2 ))- At the same time, the ratio of 
these power corrections to the total amplitude of the process remains within reasonable 
limits. To see this, we normalize the corrections T 2 in terms of the T LO ignoring the 
different tensor factors in (JU) and setting A^((,W 2 ) = A((,W 2 ), £>J(1 GeV 2 ) = 1, but 
keeping the factor J2 e l m definition of the function T (Q 2 , £, W 2 ). The ratio T 2 PC /Tq° 
calculated in this approximate way is shown in Fig. |3J We find that power corrections 
T 2 PC may amount to 31 - 38 % of the T LO at Q 2 = 1 GeV 2 and only to 6 - 7 % of its 
value at Q 2 = 10 GeV 2 . The precise estimate of the effects generated by the helicity-two 
component of the amplitude (jlj) requires more detailed investigation. 

The HT ambiguities (fTH|) coming from the principal value prescription, in the process 
under consideration are sizeable only for small values of the momentum transfer 1 GeV 2 < 
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Q 2 < 2.5 GeV 2 [from ±9% to ±3%]. At Q 2 = 5 GeV 2 they are already less than ±1% of 
the original result. As an example, the relevant curves for the ratio R^{Q 2 ) are shown in 
Fig.H 

5 Concluding remarks 

In this work we have estimated the power corrections to the amplitudes Ti(Q 2 , £, W 2 ) 
of the process 7*7 — > Tin, originating from the end-point regions z — > 0, 1. To this 
end, we have employed the symmetrized RC method combined with techniques of the IR 
renormalon calculus. We have obtained the Borel resummed expressions [Tj(Q 2 , W^ 2 )] res 
for the amplitudes and have removed IR renormalon divergences by means of a princi- 
pal value prescription. In the considered process the Borel transform of the amplitude 
Tq(Q 2 , (, W 2 ) contains an infinite number of the IR renormalon poles. Since each IR renor- 
malon pole M = n in the Borel transforms B q{g) [T (Q 2 , (, W 2 )]{u), B[T 2 (Q 2 ,(,W 2 )}(u) 
corresponds to the power correction ~ 1/Q 2n to the amplitudes, the expression (|23|). in 
general, contains power corrections ~ 1/Q 2n , n = 1,2, ...00. In numerical computations 
we have truncated the corresponding scries cit ti max — 50. As an important consistency 
check, we have proved that the result obtained within the symmetrized RC method in the 
asymptotic limit Q 2 —>■ 00 reproduces the standard prediction for the amplitudes. 

It is known that the principal value prescription generates higher twist uncertainties. 
We have shown that these uncertainties at Q 2 = 1 GeV 2 do not exceed ±10% of the 
original prediction and can be safely neglected in the region Q 2 > 5 GeV . 

Our numerical calculations have demonstrated in an admittedly method dependent 
way that the power corrections coming from the analysis of end-point regions may be es- 
sential in the region of photon virtualities Q 2 ~ a few GeV 2 . Therefore, the phenomeno- 
logical analysis of the process 7*7 — * tctt in the presently experimentally-accessible energy 
regimes should include them. 
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